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ABSTRACT 


This  report  summarizes  work  carried  out  during  1967  under  Contract 
AF49 (638)-166S  on  detonation  characteristics  of  low-density  granular 
materials.  The  purpose  of  the  research  under  this  and  preceding  contracts 
is  to  clarify  the  processes  during  initiation  to  detonation  in  propellants 
and  explosives.  Toward  this  end  studies  have  been  made  of  (a)  the 
mechanisms  of  initiation  and  propagation  in  low-density  NH^CIO^,  (b)  the 
properties  of  Initiating  shocks,  (c)  the  loci  of  shocked  states  achieved 
in  low-density  granular  aggregates,  (d)  hot  spots  in  shocked  granular 
NH^ClOj  and  KC1 ,  and  (e)  friction  and  jetting  as  possible  causes  of  hot 
spots . 

The  main  conclusions  are  that  shocks  greater  than  about  5  kbar  are 
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reactive  in  granular  NH^CIO^  at  a  density  of  1.0  g/cm  ,  that  the  reaction 
occurs  from  the  time  of  shock  entry  and  therefore  is  at  least  partially 
in  the  shock  front,  and  that  the  likely  sites  of  this  shock-front  reaction 
are  hot  spots  which  occur  when  the  interstices  collapse.  These  hot  spots 
have  been  photographed  at  about  40  kbar. 
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0.1  /nsec  was  the  sho.  test  usable  exposure  time.  A  typical  example  of  the 

hot  spots  produced  in  NHXIO,  with  the  40  to  50  kbar  shock  pressure 
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induced  is  shown  in  Fig.  4.  A  similar  picture  of  hot  spots  in  a  KC1  bed 
is  shown  in  Fig.  5. 


FIG.  4  IMAGE  CONVERTER  PICTURE  OF  HOT  SPOTS  IN  AN  AGGREGATE 
OF  NH^CIO^.  The  shock  pressure  wos  between  40  and  50  k i lobar s .  ”he 
atmosphere  was  methane.  The  Beckman  and  Whitley  image  converter 
was  operated  at  on  exposure  time  of  0.1  fisec  with  ASA  10,000  Polaroid  film. 

The  light  intensity  and  duration  measurement  was  made  through  the 
same  objective  lens  by  using  a  slightly  silvered  mirror  placed  ahead  of 
the  fast-shut ceri ng  image  converter  tube.  The  light  measurement  was  thus 
the  total  output  over  the  entire  crystal  bed  as  a  function  of  time  as  the 
shock  wave  traversed  the  depth  of  the  bed.  The  time  of  the  image  converter 
exposure  can  be  correlated  with  the  photomultiplier  record.  For  crystal 
beds  approximately  3  mm  deep,  the  light  lasted  about  2  /nsec  in  all  cases. 
When  the  crystals  were  NH^CIO^,  the  light  output  peak  was  two  to  three 
times  that  reached  when  the  crystals  were  KC1 .  This  observation  does  not 
in  itself  indicate  a  difference  in  reactivity  between  the  materials,  since 
the  KC1  crystals  are  noc  completely  transparent  and  might  be  expected  to 
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INTRODUCTION 


Tl. is  report  describes  work  carried  out  between  1  January  1967  and 
31  December  1967  on  Contract  At’  49(638) -1669,  "Detonation  Characteristics 
of  Low-Density  Granular  Materials.  1  The  purpose  of  the  work  is  to  develop 
a  basic  understanding  of  factors  controlling  the  shock  sensitivity  of 
propellants  and  explosive  materials. 

Shock  sensitivity  is  a  fundamental  property  of  propellants  or 
explosives  which  often  determines  their  vulnerability  to  unfavorable 
environments.  Detailed  analytical  and  experimental  studies  on  the  shock 
sensitivity  of  ammonium  perchlorate  (a  component  of  many  propellant 
mixtures  and  a  low-powered  insensitive  explosive  in  its  own  right)  have 
shown  a  number  of  initiation  characteristics  that  apparently  apply  to 
granular  explosives  generally.  Thus,  studies  on  this  material  are  of 
interest  not  only  to  illustrate  the  behavior  of  a  representative 
propellant  component  but  also  to  demonstrate  the  properties  of  high 
explosives.  As  a  high  explosive,  ammonium  perchlorate  (AP)  might  be 
expected  to  be  unusual  because  of  Its  low  power  and  insensitivity,  but 
the  iatter  characteristic  actually  is  particularly  favorable  for  shock 
initiation  studies.  The  shocks  necessary  for  initiation  of  low-density 
which  are  above  15  kbar,  are  easier  to  control  and  measure  than 
shocks  below  3  kbar,  which  are  adequate  to  initiate  PETN  with  measurable 
delays . 

Most  of  our  investigations  have  thus  far  been  carried  out  at  densities 
3 

of  1  g/ern  ,  that  Is,  on  pressings  of  high  porosity.  By  working  with  this 
porosity  it  is  possible  to  emphasize  the  effects  of  heterogeneity.  The 
advantage  of  the  combination  of  low  density  and  insensitivity  is  demon¬ 
strated  by  the  accessibility  of  low-pressure  shocked  states,  which  makes 
It  possible  to  define  a  region  o'  positive  slope  on  the  locus  of  shocked 
states  is)  the  p.v  pi  are. 

bast  year  our  research  was  directed  toward  demonstrating  that  the 

shock  front  in  a  granular  reactive  material  such  as  NH  CIO  was  quite 
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different  from  a  shock  front  in  an  inert  material  'ch  as  KC1  even  at 
shock  pressures  as  low  as  10  kbar.1  Sound  velocities  of  porous  compacts 
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of  these  materials  were  used  to  asaist  in  constructing  a  shock -velocity/ 
particle-velocity  plot,  which  demonstrated  a  pronounced  difference  between 
NH^CIO^  and  KC1  aggregates  above  a  particle  velocity  of  0.5  mm/Jtsec 
(corresponding  to  about  10  kbar)  .  The  difference  between  the  two 
materials  was  also  quite  obvious  in  the  p,v  plane.  Measurements  with 
a  light  gas  gen  were  used  to  confirm  previous  determinations  on  porous 
NH4C104  made  with  an  explosive  shock-generating  system.  It  remained  to 
determine  the  compressibility  of  the  KOI  aggregate  at  about  40  kbar  and 
to  determine  the  Hugoniot  of  NH^CIO^  at  crystal  density. 

Carrying  on  toward  the  general  objective  of  this  research,  which  is 
to  clarify  the  shock  initiation  and  transition-to-detonation  processes, 
we  have  this  year  worked  in  three  areas: 

1.  Confirmation  and  extension  of  the  loci  of  shocked  states 

for  NH  .CIO.  and  KC1 . 
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2.  Experimental  investigation  of  the  nature  of  hot  spots  in 
granular  materials. 

3.  Formulation  of  theories  of  hot-spot  formation  according  to 
the  mechanisms  of  friction  and  hydrodynamic  jetting. 

II  SHOCK  COMPRESSION  EXPERIMENTS 

In  last  year's  annual  report,  detailed  experimental  evidence  was 
presented  to  show  that  the  ahock  compression  curve  for  granular  ajgregctes 
of  AP  was  quite  different  from  that  for  aggregates  of  KC1  at  porosity 
m  =  VD/VC  =  2.  Two  points  remained  to  be  established:  The  Hugoniot  for 
NH^CIO^  at  crystal  density  needed  to  be  established  more  definitely  so 
that  t.ie  material  could  be  compared  with  AP  granular  aggregated  and  with 
other  solid  crystals.  The  state  variables  for  KC1  aggregates  at  m  -  2 
were  needed  for  pressures  in  the  neighborhood  of  40  k^ar  so  that  a  clear 
comparison  with  the  NHdClC>4  aggregate  could  be  made  in  this  pressure 
region. 

A .  Hugoniot  for  Single  Crystal  Ammonium  Perchlorate 

Previously  we  have  drawn  the  ammonium  perchlorate  crystal  Hugoniot 
to  match  the  known  curve  for  KC1 ,  using  sound  velocity  measurements  to 


retnforce  our  confidence  in  the  low  part  of  the  curve.  It  will  be 
desirable,  however,  to  have  detailed  information  on  the  Hugoniot  of 
NH^CIO^  crystals  to  establish  directly  the  magnitude  of  tne  deviation 
of  the  shock  compression  curve  for  the  aggregate. 

A  Hugoniot  compression  point  for  a  stngle  AP  crystal  has  been 
establtshed  using  the  light  gas  gun  techntque.  (The  method  was 
described  ir,  the  1366  annual  report  and  detatls  ot  the  technique  are 
given  in  Ref.  2.)  The  measurement  was  for  a  plane  shock  perpendtcular 
to  the  broad  face  of  the  crystal,  the  direction  in  which  the  lattice 
dimension  is  7.453  8.  Because  of  crystal  shape,  it  is  not  feasible  to 
measure  the  compression  in  other  dtractions  wtth  presently  obtainable 
crystals;  however,  the  measurement  made  is  sufficient  for  compartson 
of  this  Hugoniot  with  that  of  other  crystals  and  with  randomly  arranged, 
porous  AP. 

A  brass  plate  about  a  quarter-inch  thick  was  accelerated  in  the 

light  gas  gun  to  0.408  mm/ilsec  and  allowed  to  impact  on  a  brass  plate 

of  similar  thickness  in  the  target  assembly.  A  shock  at  about  0.20 

ram/jisec  (corresponding  to  roughly  68  kbar)  was  thus  formed  in  the 

* 

stationary  brass  plate.  An  ammonium  perchlorate  crystal  1.763  mm  thick 
was  mounted  on  the  rear  face  of  the  plate  with  a  quartz  gage  on  its  far 
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surface  to  record  arrival  times  and  allow  estimation  of  pressures  The 

double-fronted  wave  recorded  is  shown  in  Fig.  1.  By  using  the  wave 

velocities  corresponding  to  the  crystal  transit  times  indicated  by  the 

two  fronts  in  this  trace  it  is  possible  to  calculate  the  pressure  and 
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volume  of  the  final  shocked  state  in  the  AP  crystal .  The  resulting 
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state  point,  p  -  21  kbar  and  v  =  0.40  cm  /g,  i  s  indicated  on  the  plot  in 
Fig.  2.  It  lies  close  to  the  KC1  crystal  Hugoniot. 

The  nature  of  the  trace  in  Fig.  1  contains  qualitative  information 
about  tne  shock  in  crystalline  ammontum  perchlorate  which  has  not  been 
reported  before.  First  appeartng  on  the  trace  is  a  low-ar.iplttude  ramp  A 


The  crystals  were  generously  furnished  by  Dr.  William  R.  McBride  of 
the  Naval  Weapons  Center,  China  Lake,  California. 
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FIG.  1  QUARTZ  GAGE  TRACE  OF 
SHOCK  ‘N  AMMONIUM 
PERCHLORATE  CRYSTAL 

Final  pressure  about  21  kilobars 


FIG.  2  SHOCK  HUGONIOT  FOR 

AMMONIUM  PERCHLORATE 
CRYSTAL.  The  specific 
volume  of  the  crystal  is  vc. 


whose  cause  has  not  been  determined  and  which  can  probably  be  ignored  in 
view  of  its  small  magnitude.  The  small  negative  signal  at  the  beginning 
of  the  ramp  is  an  electronic  reaction  produced  by  the  particular  pre¬ 
amplifier  used.  The  second  feature  of  the  trace  is  a  rapidly  rising 
signal  B  traveling  at  4.3  mm/lfsec.  The  final  pressure  rise  is  traveling 
at  3.3  mm/fisec,  and  produces  the  final  compression  in  the  crystal.  This 
double-fronted  wave  in  ammonium  perchlorate  probably  arises  from  a 
polymorphic  transition.  The  height  of  the  precursor  B  indicates  that 
this  phase  change  may  take  place  as  low  as  5  kbar.  Since  the  existence 
of  the  new  phase  in  the  solid  crystal  is  somewhat  peripheral  to  shock 
initiation  of  granular  aggregates,  we  have  noc  performed  confirming 
experiments.  However,  the  possibility  of  this  phase  will  be  of  interest 
to  those  concerned  with  the  effects  of  shock  on  solid  ammonium  perchlorate, 
and  it  may  be  an  important  consideration  in  the  shock  vulnerability  of 
AP-containing  propellants. 
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B.  Shocked  Porous  KC1 


In  last  year's  annual  report,  measurements  on  aggregates  of  K'Cl  at 

porosity  m  =  v^/v^  «  2  were  reported  in  the  pressure  range  below  20  kbar. 

Measurements  were  made  with  explosive  shock-generating  systems  and  with 

the  light  gas  gun.  It  was  apparent  that  the  shock-compacted  state  of 

these  aggregates  lay  close  to  the  KC1  crystal  Hugoniot.  Although  it  is 

thermodynamically  impossible  for  them  to  lie  exactly  on  the  crystal  curve 

their  position  is  very  close  to  that  called  for  by  the  compressible 
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locking -solid  model.  Linde  and  Schmidt  report  such  dynamic  results  for 
porous  aluminum  samples  (m  =  1.3  to  m  =  2.5).  Furthermore,  samples 
recovered  after  shocking  to  pressure  in  the  range  from  5  to  60  kbar 
exhibited  densities  close  to  that  of  solid  aluminum.  Dynamic  measure¬ 
ments  also  showed  that  the  compacted  volumes  of  porou3  graphite  above  a 

few  kilobars  were  essentially  those  of  the  solid  material  at  the  same 
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pressure.  On  the  other  hand,  Russian  data  on  porous  KC1  at  higher 
pressure  indicate  that  the  shocked  states  of  the  aggregates  occupy 
noticeably  greater  volumes  than  those  of  the  solid  crystal.  The  two 
sets  of  results  are  not  in  direct  contradiction,  since  the  pressures  are 
quite  different.  But  it  is  important  for  our  purposes  to  know  to  what 
pressure  the  complete  collapse  of  the  KC1  aggregate  persists  and  to  be 
able  to  define  the  nonreactive  KC1  curve  ovev  the  region  in  which  the 
specific  volume  of  shocked  NlljClO^  aggregate  is  larger  than  crystal 
density . 

Two  experiments  were  conducted  with  an  explosive  shock -generating 
system  using  the  techniques  described  at  the  Fourth  Detonation  Symposium. 
A  pressure  in  the  neighborhood  of  40  kbar  was  obtained  with  an  explosive 
shock  system  consisting  of  19  mm  of  magnesium  and  9  mm  of  Plexiglas  on  a 

P-80  plane  wave  generator.  The  aggregate  beds  were  pressed  from  KC1 

3 

particles  with  diameters  of  IS  /j  to  a  density  of  ]  .0  g/cm  (m  2)  , 

Seven  samples  were  shucked  in  the  first  experiment  and  five  in  the  second 
The  results  of  the  two  experiments  were  essentially  in  agreement,  as 
shown  by  the  open  squares  near  40  kbar  in  Fig.  3. 

The  conclusion  from  these  experiments  is  that  KC1  aggregates  at 
porosity  m  =  2  collapse  essentially  to  the  volume  of  compressed  crystal 
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FIG.  3  SHOCK  COMPRESSION  OF  POTASSIUM  CHLORIDE  AGGREGATES 

AT  POROSITY  m  =  2.  New  data  are  the  squares  at  the  poinrs  of  arrows 
originating  at  initial  specific  volume  of  the  aggregates  vQ. 


KC1  at  the  same  pressure.  We  must  accept  the  data  of  Kormer  et_  aJL 
indicating  that  at  higher  pressures  the  shocked  states  are  noticeably 
different  from  those  attained  by  the  shocked  crystal.  Apparently  this 
difference  becomes  noticeable  somewhere  between  43  and  89  kbar. 

HI  HOT  SPOTS  IN  SHOCKED  GRANULAR  MATERIALS —EXPERIMENTAL  WORK 

In  our  work  on  the  shocked  states  of  AP  aggregates  we  have  shown 
that  the  locus  of  shocked  states  Is  different  tn  kind  from  the  locus 
for  KC1  aggregates.  These  loci  have  been  constructed  on  the  basis  of 
the  entering  shocks  m  the  aggregates.  Therefore  we  have  concluded  that 
a  certain  amount  of  chemical  reaction  occurs  at  the  front  of  the  shock 
in  AP  aggregates,  and  that  this  accounts  for  the  different  shocks  state 
loci.  If  such  reaction  is  to  occur  at  the  shock  front  it  will  occur  at 

the  hottest  regions  there.  Hot  spots  have  been  observed  in  detonating 
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materials,  in  initiating  shocks,  and  in  shocks  in  inert  materials.  1 

We  have  therefore  undertaken  a  study  of  the  characteristics  of  the  hot 

spots,  to  establish  the  connection  between  these  centers  of  early  reaction 

in  the  shock  wave  and  the  displaced  locus  of  shocked  states  in  granular 

ammonium  perchlorate. 

We  have  successfully  photographed  hot  spots  in  shocked  granular 
materials  with  a  Beckman  and  Whitley  image-converter  camera  at  very 
short  exposure  times.  Beds  of  crystals  1  to  2  mm  in  size,  both  KOI  and 
NH^CIO^,  were  shocked  by  a  calibrated  svstem  consisting  of  an  explosive 
pellet  and  a  Lucite  attenuator.  The  Lucite  was  thoroughly  blackened  to 
eliminate  light  from  the  pellet  and  the  bed  was  flooded  with  methane  to 
avoid  gas  shock  lumir  _ity.  Relative  intensity  and  duration  of  the  light 
generated  by  the  crystals  were  recorded  with  a  photomultiplier  tube. 

This  exploratory  phase  of  the  study  has  covered  two  pressure  ranges, 
the  higher  being  somewhat  above  100  kbar  and  the  lower  being  approxi¬ 
mately  40  to  50  kbar  in  the  crystal  bed.  Precise  statement  of  pressure 
in  this  ^arse-grained,  heterogeneous  material  is  not  possible  using 
continuum  hydrodynamic  methods.  Successful  ohotographs  of  the  hot  spots 
were  made  at  the  high  pressure  range  with  exposures  as  short  as  0.05  iisec. 
At  the  low  pressure  range,  with  less  light  produced  in  the  shock  process, 
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0.1  fiscc  was  the  shoi  test  usable  exposure  time.  A  typical  example  of  the 
hot  spots  produced  in  NH^CIO^  with  the  40  to  50  kbar  shock  pressure 
induced  is  shown  in  Fig.  4.  A  similar  picture  of  hot  spots  in  a  KC1  bed 
is  shown  in  Fig.  5. 


FIG.  4  IMAGE  CONVERTER  PICTURE  OF  HOT  SPOTS  IN  AN  AGGREGATE 
OF  NH^CIOj.  The  shack  pressure  was  between  40  and  50  kilobars.  ~he 
atmosphere  was  methane.  The  Beckman  and  Whitley  image  converter 
was  operated  at  an  exposure  time  af  0,1  psec  with  ASA  10,000  Polaroid  film. 

The  light  intensity  and  duration  measurement  was  made  through  the 
same  objective  lens  by  using  a  slightly  silvered  mirror  placed  ahead  of 
the  fast-shut cering  image  converter  tube.  The  light  measurement  was  thus 
the  total  output  over  the  entire  crystal  bed  as  a  function  of  time  as  the 
shock  wave  traversed  the  depth  of  the  bed.  The  time  of  the  image  converter 
exposure  can  be  correlated  with  the  photomultiplier  record.  For  crystal 
beds  approximately  3  mm  deep,  the  light  lasted  about  2  /jsec  in  all  cases. 
When  the  crystals  were  NH^CIO^,  the  light  output  peak  was  two  to  three 
times  that  reached  when  the  crystals  were  KC1 .  This  observation  does  not 
in  itself  indicate  a  difference  in  reactivity  between  the  materials,  sines 
the  KC1  crystals  are  noi:  completely  transparent  and  might  be  expected  to 
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attenuate  and  scatter  light  more  than  the  NH  CIO  crystals.  However. 

4  4 

since  we  have  previously  observed  a  dramatic  difference  in  shock  com¬ 
pression  behavior  for  the  two  materials  when  they  are  in  finely  divided 

porous  form,  th  increased  light  output  in  shocked  NH  CIO  crystals  seems 

4  4 

most  likely  to  be  .  nother  evidence  of  early  reaction. 


FIG.  5  IMAGE  CONVERTER  PICTURE  OF  HOT  SPOTS  IN  POTASSIUM 
CHLORIDE  AGGREGATE.  Sume  conditions  os  for  Fig.  4. 

IV  HOT  SPOTS  BY  FRICTION  OR  JETTING 

We  conclude  lrom  our  work  on  the  loci  of  states  achieved  in  NH.C10 

4  4 

and  KC1  aggregates  that  reaction  occurs  in  the  front  of  initiating  shocks. 

Postulation  of  hot  spots  is  required  to  reach  reactive  temperatures. 

Very  hot  localized  regions  have  been  observed  in  detonation  fronts  in 
,  ,  6,7 

granular  explosives.  In  the  preceding  section  of  this  report  image 
converter  pictures  of  hot  localized  regions  were  shown  in  shocked  beds 
of  NH.C10,  and  KC1  granules  when  the  shock  was  around  40  kbar.  Presumably, 
cooler  hot  spots  also  exist  at  lower  shock  pressures,  and  it  is  hoped 
this  can  be  demonstrated  by  photomultiplier  measurements  over  a  lower 
range  of  pressures. 
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In  detonations,  where  the  products  can  reasonably  be  presumed  to  act 

6 

like  gases,  stagnation  is  a  probable  mechanism  for  hot  spots.  However, 

in  initiating  shocks  it  is  not  clear  that  collapse  of  the  interstices 

produces  gas  until  after  local  reaction  occurs.  It  seems  reasonable, 

therefore,  that  a  process  involving  the  solid  material  should  be 

considered.  The  following  account  summarizes  a  treatment  by  William 

10 

Band  of  local  heating  in  shocked  granular  materials  . 

Consider  a  shock  front  arriving  at  the  free  surface  of  a  cavity  in 
a  granule  (or  between  two  granules).  The  shock  reflects  at  the  free 
surface,  but  continues  on  either  side  of  the  cavity. 

The  diagram  in  Fig.  6  gives  an  idealized  geometry.  If  u  is  the 
particle  velocity  behind  the  shock,  the  velocity  of  the  free  suiface  in 
the  state  leit  behind  by  the  reflected  rarefaction  is  about  2u.  A 
relative  motion  therefore  exists  between  the  material  A,  as  it  proceeds 
to  fill  the  cavity,  and  the  material  B  at  the  sides.  If  a  pressure  is 
maintained  by  the  shock  across  this  interface,  friction  can  dissipate 
some  of  the  mechanical  energy  into  heat. 


FIG.  6  IDEAIIZEO  DIAGRAM  OF  THE  ACTION  OF  A  SHOCK 
ON  A  CAVITY  IN  A  CRYSTAL  OR  IN  AN  INTERSTICE 

Temperatures  are  calculated  as  the  result  af  friction  between 
material  from  the  tap  of  the  cavity  moving  at  twice  the  particle 
velocity  2u  and  material  above  the  side  walls  moving  at 
velocity  u. 
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In  many  materials,  solid  rigidity  does  not  vanisi.  immediately  but 
may  persist  for  a  time  from  10  nsec  to  as  long  as  1  p sec.  After  the 
onset  of  the  shock  pressure  there  will  be  a  very  short  time  T  during 
which  it  would  be  legitimate  to  think  of  the  friction  as  sliding  friction 
between  two  solid  surfaces.  During  this  short  time,  shear  stresses 
comparable  with  the  pressure  p  will  do  work  by  friction  over  the 
relative  displacement  of  the  two  surfaces. 

Actually,  the  shock  velocity  is  usually  considerably  greater  than 
2u,  so  the  high-pressure  slippage  is  limited  to  the  region  above  the 
cavity.  We  approximate  the  area  of  effective  contact  A  =  JfsUT  where  s 
is  the  linear  dimension  of  the  cavity  and  U  is  the  shock  velocity. 

The  distances  of  relative  motion  between  the  surfaces  is  uT.  The  work 
done  by  friction  is  thus 

2  2 

A  W  =  0pAu7  =  iSp  ST  /p 

where  0  is  the  coefficient  of  friction,  the  volume  heated  is  AX,  X 
is  the  thickness  of  the  heated  layer,  and  we  assume  this  is  determined 
by  thermal  diffusion  during  the  time  T  Lambda  is  then  given  by 
X  =  ^TTKT/pC  where  K  is  the  thermal  conductivity  and  C  is  the  specific 
heat.  The  maximum  temperature  rise  from  this  source  of  heat  is  therefore 

A  T  -  f»puT'pcX  =  OpuyJrS  ^ITKfZ  (1) 

After  a  time  t  >  T,  the  temperature  drops  by  diffusion: 


Some  modification  of  this  statement  is  to  be  expected  from  the  geometry 
of  the  cavity  surface,  but  as  an  order  of  magnitude  calculation,  (1)  and 
(la)  should  give  a  good  idea  of  expected  heating.  There  is  nothing  in 
this  result  to  account  for  the  significant  effects  of  cavity  size  that 
are  observed  experimentally. 

As  the  shock  envelops  the  cavity,  the  sideways-propagating  relief 
waves  develop  and  the  cavity  collapses  from  the  sides.  Material  from  the 
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boundary  layer  between  A  and  Bis  thus  caused  to  jet  into  the  cavity  and 
may  strike  the  floor  of  the  cavity  before  it  has  collapsed.  If  this 
impact  occurs  on  previously  unshocked  material,  the  same  sudden  shear 
stresses  can  be  initiated  at  the  boundaries  of  the  area  struck  by  the 
hypervelocity  jet.  Another  increase  in  temperature  therefore  occurs 
here,  Riven  by  the  same  type  of  relation: 

At2  =  fjp'u  yfr/yfmpc  (2) 

•  / 

where  p  and  u  refer  to  the  shock  induced  in  the  floor  of  the  cavity 
by  the  Jet  impact . 

*  * 

The  ratio  p  u  /pu  depends  on  the  cavity  size,  and  we  shall  develop 
a  simple  model  for  which  this  is  true.  We  therefore  accept  the  ratio 

jS  =  0p  'u  V pu 

as  an  effective  conversion  factor  and  write 


Results  are  listed  in  Tables  I  through  III  for  PETN,  KC1 ,  and  NaCI .  In 
using  (2a)  we  are  plagued  with  several  unknown  parameters,  the  >-orst  of 
which  is  the  "coefficient  of  friction,”  or  the  ratio  of  the  shear  stress 

to  the  pressure  during  initial  slipping.  If  0  is  assumed  as  small  as 

*  *  —  8 
O.Oi,  p  u  /pu  no  greater  than  5,  and  T  =  10  sec,  the  temperature  in 

PETN  will  come  up  to  A  =  900°. 

The  next  problem  is  to  consider  the  relation  between  /3  and  the 

cavity  size. 

In  these  calculations  we  nave  assumed  that: 

I,  The  rarefaction  front  travels  a  distance  OC  =  ut  (Fig.  6) 
during  the  operation  of  the  friction  before  it  arrives  at  the 
back  of  the  granule: 

s  >  uT 

e 


12 


TABLE  I 


LOCAL  HEATING  IN  GRANULAR  PETN3 


P/109  2 
(dynes/cm  ) 

u/104 

(cm/sec) 

(s  /s  )* 

8  C 

A  T/103/9 
Deg  .C 

A  TjS 
Deg  .C 

1 

0.4 

2.74 

0.32 

12.9 

2 

0.75 

2.06 

1.21 

45.0 

3 

1.1 

1.69 

2.66 

89.0 

4 

1.45 

1.53 

4 .67 

170.0 

5 

1.8 

1.43 

7.25 

260.0 

6 

2.1 

1.38 

10.1 

353.0 

7 

2.35 

1.35 

13.3 

442.0 

8 

2.6 

1.33 

16.7 

541.0 

» 

2.85 

1.31 

20.7 

650.0 

10 

3.1 

1.29 

25.0 

769.0 

12 

3  .5 

1.27 

33.8 

980.0 

14 

3.9 

1.26 

44.0 

1220.0 

16 

4.15 

1  .26 

53.5 

1380.0 

18 

4.4 

1 .26 

63.8 

1550.0 

20 

4  65 

1 .26 

74.9 

1730.0 

-3 

3  =  1 . 8  g  cm  , 

K  =  2.5xl04  g 

-3  • 

cm  sec  deg 

-1.  C  =  1 

7 

.1x10  erg 

T  =  10  8sec,  1/ =  8.05x10”',  X  =  -^77TK/pC  =  6.3x10 


b.  This  column  gives  the  temperatures  expected  if  we  have  only  the 
heat  obtained  by  ’stagnation"  of  the  jet  at  speeds  about  4u  — 
A  T  =  i(4u)2/C. 


deg 
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TABLE  II 


LOCAL  HEATING  IN  GRANULAR  KCl* 


p/109  2 
(dynes/cm  ) 

4 

u/10 

(cm/sec) 

(s  /s  ) 

e  c 

A  T/1O30 
Deg  .C 

3° 

A  Tj/10 

Deg  .C 

10 

1.5 

2.5 

2 . 83 

0.264 

20 

3.0 

1.75 

11.34 

1 .06 

30 

5.8 

1.33 

32.9 

3.95 

40 

6.8 

1 .31 

51.4 

5.43 

eo 

8.8 

1 .26 

99.8 

9.10 

80 

10.7 

1 .24 

162.0 

13.45 

100 

12.9 

1.20 

244.0 

19.55 

3  5 

p  =  J  984  g/cm  ,  K  =  6.7x10  g  cm  sec 

3  -1 

deg  ,  C  =  0 

,68xl07erg  deg 

-  o 

T  =  10  sec, 

1/yjmfC  =  l 

.89x10  ,  X  = 

yimc/pc  =  1.25x10 

This  column 

gives  the  rise  in  temperature  expected 

if  only  the 

heat  obtainable  from  the 

cavity-jet  at 

speed  4u  is 

included . 
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TABLE  III 


LOCAL  KEATING  IN  GRANULAR  NaCI® 


p/109  2 

(dyne  s/cm  ) 

u/104 

(cm/sec) 

(s  /s  ) 

g  c 

*  At/ip3(S 

Deg  .C 

A  T, 

Deg  .C 

2 

0.3 

6.2 

0.126 

8.4 

4 

0.55 

4.1 

0.462 

28.0 

6 

0.80 

3.2 

1.080 

60.0 

8 

1.05 

2.7 

1.760 

103.0 

10 

1  .13 

2.8 

2  .370 

120.0 

12 

1.55 

2.2 

3  .910 

224.0 

14 

1.80 

2.0 

5.290 

303.0 

16 

2.00 

1  .9 

6.820 

374.0 

18 

2.25 

1.8 

8.510 

473.0 

20 

'1.45 

1.8 

10.290 

560.0 

30 

3.60 

1.5 

22.700 

1210.0 

p  =  2.165  g/ 

3  5 

cm  ,  K  =  4.0x10 

-3.  -1  „ 

g  cm  sec  deg  ,  C  u 

0.84xl07erg  g  ^deg 

T  =  10  *sec, 

1/^TOpC  =  2.1xl0"7,  X 

=  yJrrK/pC  =  2. 

63  xl0~5 

This  column 

gives  the  temperature  rise  expected  if 

only  the  heat 

obtained  by 

"stagnation"  of 

the  jet 

at  speed  about 

4u  is  included. 
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2.  The  reflected  shock  from  the  bottom  of  the  cavity  does  not 
limit  the  depth  of  the  low-pressure,  high  particle-velocity 
material  to  be  less  than  the  distance  ut: 


s  (1  +  U/2u)  >  UT 
c 

If  either  or  both  of  these  conditions  are  violated,  the  correct  formulas 
become  the  fol. owing: 


If  s^  <  uT,  friction  operates  for  only  the  time  T 
the  work  done  is  only 


s  /u  and 
g 


A W  =  0pAs 
1  g 

while  the  temperature  increase  is 

At’  =  flpsg/X’pc  (ib) 

where 

X'=^JTKT  '/pC 

or  finally 

A  t'  =  0p  us^/^TKpC  (lc) 

If  s  <  uT/(l  +  U/2u)  <  s  ,  heating  is  reduced  by  the  factor  s  /s  and 
c  g  eg 

the  correct  relation  is 


Similar  modifications  would  be  expected  in  the  expressions  for  At^. 

If  s  >  ur  while  s  <  uT(l  +  U/2u) ,  we  must  return  to  (la)  and  reduce 
g  c 

that  by  the  factor  s^d  +  U/2u)/uT: 

ATj”  --  0psc(l  +  U/2u)/^f  TTKpe  (le) 

These  results  suggest  that  experiments  should  be  designed  to  find  the 

value  of  T.  For  a  fixed  porosity  ratio  (s  /s  fixed)  find  the  maximum 

eg 

grain  size  beyond  which  the  temperature  no  longer  appears  to  increase 

with  s  ; 
g 
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T  =  s  /u 

gmax 

Then  lor  a  fixed  s  <  s  look  lor  a  maximum  porosity  or  cavity  size 

g  gmax 

beyond  which  the  temperature  no  longer  increases  with  s^; 

T  =  (1  +  U/2u)s  /u 
cmax 

If  the  cavity  size  s^  is  less  than  ut ,  there  is  not  sufficient 
time  for  the  material  to  lose  its  solid  rigidity,  but  if  s  >  uT,  then 
the  solid  behaves  hydrodynamically ,  and  material  can  jet  into  the  cavity. 

Consider  a  cylindrical  cavity  in  which  the  shock  has  been  reflected 
from  the  top  (Fig.  7).  This  surface  is  moving  down  at  speed  2u,  and 
the  sides  are  moving  ir.  at  speed  a.  The  angle  of  the  interface  is 
0  =  tan  (more  exactly  it  is  greater  than  this,  because  the  sides  are 

inclined  to  the  axis  of  the  cylinder  because  of  the  finite  speed  of  the 
shock).  To  reduce  the  point  of  contact,  0,  to  rest,  we  impress  a 
velocity  -yiu  along  the  line  of  contact.  The  flow  picture  is  shown  in 
Fig.  M.  Bernoulli's  equation  for  uniform  pressure  (in  the  cavityl 
evidently  requires  the  jet  to  be  a  composite  of  two  parts  having  different 
speeds  v  and  2u.  We  write  Of  for  the  fraction  of  the  material  arriving 
from  the  top  of  the  cavity  that  goes  into  the  jet,  and  j3  for  the  corres¬ 
ponding  fraction  irom  the  sides.  Conservation  of  momentum  leads  to  the 
results 

a  =  id  -  l/\fs)  0  =  i(l  -  2/y[T) 

=  0.277,  =  0.051. 

We  now  have  to  picture  a  cone  of  materials  (shown  in  Fig.  9)  with  a 
vertex  at  rest  in  the  laboratory  coordinates  as  the  cavity  collapses. 

The  inside  of  the  cone  is  flowing  at  the  speed  u  +  ^5u,  and  the  outer 
part  of  the  cone  at  the  speed  2u  »  ^u.  T  >  conserve  momentum  it  turns 
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FIG.  7  CONDITIONS  AT  THE  CORNERS  FIG.  8  JETTING  AT  THE 

OF  A  CYLINDRICAL  CAVITY  CORNERS  OF  A 

UNDER  THE  ACTION  OF  A  CYLINDRICAL  CAVITY 

SHOCK  WITH  PARTICLE 
VELOCITY  u  NORMAL  TO  THE 
TOP  OF  THE  CAVITY 


•l 

li-lWOI  •Z-ZO* 


FIG,  9  INTERACTION  OF  JETS  IN 
A  CYLINDRICAL  CAVITY 
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out  that,  for  the  present  geometry,  the  internal  jet  forms  entirely  from 
the  inner  layer,  while  the  down-jet  is  composite:  all  the  outer  layer 
goes  down  and  a  core  from  the  inside  of  the  cor.e  also  goes  down. 
Conservation  of  momentum  leads  to  the  result  that  the  fraction  of  the 
inner  layer  jetting  up  is 

Ctj  =  (1  -  sin6)  (i  +  ^u2/20u1) 

where  u^  =  u(l  +^5)  ,  u^  =  u  (2  +^5)  .  Numerically  this  turns  out  to 

be  a  =  0.0675 . 

J 

The  net  result  is,  then,  that  a  fraction  Q  x  ^  =  0.019  of  the 
material  in  the  top  of  the  cavity  jets  back  up,  while  (1  -  Qf)Qi  =  23 
of  the  material  from  the  top  and  0.054  of  the  material  in  the  sides 
jets  down  in  a  composite  motion — the  material  from  the  sides  having  a 
speed  4.24  u  and  that  from  the  top  (in  the  core  of  the  jet)  only  a 
speed  3.24  u.  Assuming  that  this  tidy,  streamline  picture  can  hold 
steady,  we  can  imagine  the  head  of  the  jet  punching  into  the  bottom  of 
the  ca"ity  and  causing  a  ring-shaped  hot  spot  by  the  same  mechanism 
already  postulated  for  the  top  of  the  cavity.  Both  the  pressure  and 
the  particle  velocity  induced  ir  the  solid  would  in  first  approximation 
be  about  double  those  in  the  original  shock,  so  a  rough  estimate  of 
p  uVpu  is  a  factor  of  4.  It  seems  more  realistic  to  suppose  that  the 
composite  nature  of  the  flov  both  in  the  cone  and  in  the  jets  would 
induce  turbulent  degeneration  of  energy,  and  that  the  jet  would  there¬ 
fore  give  a  much  smaller  punch  to  the  base  of  the  cavity  than  estimated 
above . 

In  spite  of  the  crudeness  of  the  model  and  the  .-’ssumptions  that  are 
necessary,  the  relationships  between  grain  size  and  cavity  size  might 
give  an  insight  into  the  consistency  of  the  model.  The  ratio 

s  /s  =  1  +  U/2u 
gmax  cmax 

should  be  amenable  to  experimental  check.  Even  though  the  absolute 
values  of  the  area  of  friction  may  be  wrong,  the  relative  effects  of 
grain  and  pore  size  should  be  as  indicated. 
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VI  PUBLICATIONS 

During  1967  the  following  research  note  was  published: 

Reese,  B.O.,  J.H.  Blackburn,  L.B.  Seely,  and  Marjorie  W.  Evans, 
"Longitudinal  Sound  Velocities  of  Granular  Compacts  of  Ammonium  Peru:  -orate 
and  Potassium  Chloride,"  Combustion  and  Flame  11,  262-264  (1967), 
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This  report  summarizes  work  carried  out  during  1967  under  Contract  AF49(638)- 
1669  on  detonation  characteristics  of  low-density  granular  materials.  The  purpose 
of  the  research  under  this  and  preceding  contracts  is  to  clarify  the  processes 
during  initiation  to  detonation  in  propellants  and  explosives.  Toward  this  end 
studies  have  been  made  of  (a)  the  mechanisms  of  initiation  and  propagation  in  low- 
density  NH4CIO4,  (b)  the  properties  of  initiating  shocks,  (c)  the  loci  of  shocked 
states  achieved  in  low-density  granular  aggregates,  (d)  hot  spots  in  shocked 
granular  NH4CIO4  and  KC1 ,  and  (2)  friction  and  Jetting  as  possible  causes  of  hot 
spots . 

The  main  conclusions  are  that  shocks  greater  than  about  5  kbar  are  reactive 
in  granular  NH4CIO4  at  a  density  of  1.0  g/cin  ,  that  the  reaction  occurs  from  the 
time  of  shock  entry  and  therefore  is  at  least  partially  in  the  shock  front,  and 
that  the  likely  sites  of  this  shock-front  reaction  are  hot  spots  which  occur  when 
the  interstices  collapse.  These  hot  spots  have  been  photographed  at  about  40  kbar. 
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